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ABSTRACT 

If sterile neutrinos exist and form halos in galactic centers, they can give rise to 
observational consequences. In particular, the sterile neutrinos decay radiatively 
and heat up the gas in the protogalaxy to achieve hydrostatic equilibrium, and 
they provide the mass to form supermassive blackholes. A natural correlation 
between the blackhole mass and velocity dispersion thus arises \og{MBHj/MQ) = 
alog(cr/200km s"^) + (3 with a ^ 4 and (3^8. 

Subject headings: Galaxies, Sterile Neutrinos, galactic center, supermassive black- 
holes, velocity dispersion 

1. Introduction 

Understanding the nature of dark matter remains a fundamental pro blem in astrophysics 



and cosmology. Sin ce the discovery of neutrinos' non-zero rest mass (iFukuda et al.l Il998 



Bilenky et al.lll998l ). the possibility that neutrinos contribute to cosmological dark matter 
has become a hot topic again. In particular, the sterile neutrinos are a class of candidate 
dark matter particles with no standard model interaction. Although the recent MiniBooNE 
result challenges the LSND result that suggests the existence of eV scale sterile neutrinos 



(lAguilar-Arevalo et al.l 120071 ). more massive sterile neutrinos (eg. keV) may still exist. The 
fact that active neutrinos have rest mass implies that right-handed neutrinos should exist 
which may indeed be massive sterile neutrinos with rest mass m^. 



Re cently, it was proposed that a degenerate sterile n eutrino halo exists in the galactic 
center ( VioUier et al.lll993l : iMunyaneza and VioUierl |2002| ). Later, this model is ruled out 
with the availability of more precise data, particularly those on orbits of stars such as S2 



- 2 - 



near the Milky Way center (jSchodel et al.ll2002l ). Nevertheless, the existence of such a sterile 
neutrino halo at the centers of protogalaxies may still be possible. Sterile neutrinos decay 
into lighter neutrinos and photons, and provide mass to fuel the growth of supermassive 
blackholes. However, since most of these sterile neutrinos may have either decayed or fallen 
into the supermassive blackhole, the total mass of the sterile neutrino halo at the galactic 
center becomes very small (-C lO^M©) at present. Moreover, the exi stence of decaying ste rile 
neutrinos may help to solve the co oling flow problem in clust ers (jChan and Chul 120071 ) as 
well as reionization in the universe (IHansen and Haimaru |2004J ) . Therefore it is worthwhile 
to discuss the consequences of the existence of massive sterile neutrinos at the centers of 
protogalaxies, which decay into light neutrinos and photons. 

On the other hand, recent observations have led to some tight relations between the 
central black hole masses Mbhj and velocity dispersions a in the bulges of galaxies. These 
relations can be summa rized as \o^{Mbhj I M q^) = a log('(T/2 00 km s~^) + /3, where a is 



found to be 3.75 ± 0.3 (IGebhardt et al. 
by two different groups. 
a = 4.02 ± 0.32 and (3 



Tremaine et al. 



2000h and 4.8 ± 0.5 JPerrarese and MerrittI boooh 



(I2OO2I ) reanalysed both sets of data and obtained 
These results indicate that blackhole formation 



8.13 ± 0.06 

may be related to galaxy formation, which challenges existing galaxy formation theories 
dAdams et aUboOll ). 



This relation has been d e rived in recent theoretical model s ( Adams et al.ll200ll : iMacMillan and Henriks* 
2OO2I : [Robertson et ahlbood : iMurrav et allbood : iKineJ boosl . I2OO3I). We assume that a de- 
generate sterile neutrino halo exists in the center of a protogalaxy. There are two different 
deca y modes for sterile neutrinos u,^. The rn ajor decaying channel is Us Su with decay 
rate (IBarger et al.lll995l : iBoyarsky et al.ll2008l ) 



3u 



3847r3 



sin^ 2eml = 1.77 x lO'^^ sin^ 26 ( 

VlkeV/ 



' s , 



where Gp and 6 are the Fermi constant and mixing angle of sterile neutrino wit h active neu- 
trino s respectively. The mi nor decaying channel is z/^ — > z/ + 7 with decay rate tearger et al. 
I995I : ISovarskv et al.l[2008h 



9aGp .20/1 5 
— sm^ 29ml 



1.38 X 10-22 sin^ 2^ 



-V 



1 



31/1 



(2) 



10247r4""^ ^" " "'"""VlkeV/ 128 

where a is the fine structure constant. It is quite difficult to detect the active neutrinos 
produced in the major decaying channel. Therefore, we focus on observational consequence 
of the radiative decay of the sterile neutrinos with rest mass rris > 10 keV. They emit high 
energy photons (~ ms/2) which heat up the surrounding gas so that hydrostatic equilibrium 
of the latter is maintained. The sterile neutrino halo also provides mass to form a super- 
massive blackhole from a small seed blackhole. Without any further assumption, a ^ A is 
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consistent with the range of the decay rate obtained by observational data from coohng flow 
clusters. In this article, we first give a brief review on three popular analytic models that 
explain the Mbhj ~ c" relation. Then we will give a detailed description of our model and 
compare it with other existing models. 



2. Review on models of the Mbhj — relation 
2.1. Super Eddington Accretion model 



Kingl (120031 ) presented a model to explain the Mbhj — ct relation. He as sumed that the 
gas density profile of a protogalaxy is isothermal (p ~ r~^) (jKing||2003l . l2005l ) . Therefore the 
gas mass inside radius R is: 



M{R) 



R 

4:71 j pr^dr 
'o 



2fyR 

G ' 



(3) 



where fg ~ 0.16 is the cosmological ratio of baryon to total mass, assumed to be the same in 
a galaxy, and the Virial Theorem is used. Consider a super-Eddington accretion onto a seed 
blackhole. The accretion feedback produces a momentum-driven superbubble that sweeps 
ambient gas into a thin shell which expands to the galaxy. The equation of motion is 



d p. ^ GM{R)[MBH{t) + M{R)] 



^edd 



(4) 



where L^dd = 4:nGMBH{t)c/ n, with k the opacity and MBnit) the mass of the central 
blackhole at time t. Integrating twice and assuming R ^ GMbhj/'^'^, one gets: 



R^ = R'o + 2RM - 1 



MBH{t) 

M, 



(5) 



where Rq = R ai R = Rq, with Rq some large radius GMbej/c^), and M^ = fgHa^/nG"^ 



Therefore, the maximum radius R„ 

R^ 



"'max 



Rl 



is given by 



1 + 



Rl 



2a\l-MBH{t)/M, 



(6) 



When MBnit) approaches M^., -Rmax becomes very large such that the cooling of the shocked 
wind is inefficient as the cooling time t cooling oc R? and the accretion is stopped as the shell 



can escape the galaxy entirely by gas pres sure (King 



mass supply (such as in a merger), we get (lKingll2005l ) 



20031 ). Therefore, given an adequate 



BHJ 



TlG^ 



(7) 
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Here, the proportionality constant fgn/nG'^ lies within the observational constraints. To 
summarize, the Mbhj — o' relation is obtained with three important assumptions: (1) isother- 
mal gas density distribution throughout the galaxy formation, (2) super Eddington accretion, 
and (3) an adequate mass supply. 



2.2. Self-similar model 



MacMillan and HenriksenI (120021 ) obtained a relation between Mbhj and a by assuming 
that the density and velocity distributions of matter are self-similar. They assumed that the 
galaxy is formed by the extended collapse of a halo composed of collisionless matter. The 
central blackhole is grown proportionally to the ha lo as matter continues to fall in. The 
relation is given by (IMacMillan and HenriksenI |2002| ) 

'M/a-2 



logM, 



BHJ 



oc 



log a. 



(8) 



5 /a -I 

where 5 and a are scales in space and time respectively, and their ratio is related to the 
power-l aw index of the initial density perturbation e in the spherical infall model of halo 
growth flHenriksen and Widrowlll999l ): 



a 6 \ e 



(9) 



The power-law index e = (n + 3)/2, where n is the index of the primordial matter power 
spectrum P{k) oc k^. If n = —2, Eq. (8) agrees with the observation Mbhj oc a'^. This 
model involves a relation Eq. (9) which is quite model dependent. 



2.3. Ballistic model 



Adams et al.l (120011 ) assume the dark matter and baryons to be unsegregated and the 
isothermal initial mass density distribution {Mt{r) oc r). The specific orbital energy is 
conserved when the particles fall into the small seed blackhole: 



F-\?^^' no^ 



where Vr and j are the radial velocity and angular momentum per unit mass. When the 
particles fall into the equatorial plane, their pericenters are 

^ 2GMt(r) 2(t8 ' ^ ' 
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where Q is the average angular speed of the particles. In the early stage, all the particles 
fall into the blackhole until the blackhole mass reaches a critical point that corresponds to 
p = 4:Rs, where MrhU) = Mt(Rs), with -R ., = 2GMBHif) / the Schwarzschild radius. This 
gives the final relation (lAdams et al.ll200ll ) 



BHJ 



GcQ' 



(12) 



This model is based on the assumption of the isothermal distribution of matter, and there 
is a free parameter Q, which is assumed to have the same value for all galaxies. 



3. Decaying sterile neutrino halo model 



Suppose the sterile neutrino halo dominates the mass in the proto-galactic center, most 
of the mass in the blackhole comes from t he sterile neutrino halo with radius R, which was 
formed in the very early universe t ~ (IMunyaneza and Biermanru l2005l : I Chan and Chu 
20071 ). The total mass of the degenerate sterile neutrino halo at time ti, is 



Jo 



(13) 



where ps is the mass density of the sterile neutrino halo. Assume that a seed blackhole with 
mass MsHfl of order solar mass is formed at tj,, long after the formation of the sterile neutrino 
halo. It would grow by accreting mass of the sterile neutrino halo to mass MBH{t) at time t. 
As some sterile neutrinos would be accreted by the seed blackhole, the degenerate pressure 
is decreased and more st erile neutrinos will fall into the b lackhole as their Fermi speed is less 
than their escape speed (IMunyaneza and Biermannll2005). The falling time scale at distance 
r from the blackhole in a free falling model is given by (jPhillipslll994l ) 



TT 



2y 2G[MBHfi + Ms{r,t)y 



(14) 



For Ms{t) > lO^Mo, > 10 keV and R < 0.04 pc, tff < 160 years for r < R, which 
is much shorter than the Hubble time. Therefore, we do not need any intermediate mass 
blackholes since the small seed blackhole can grow to a 10^ — IO^Mq supermassive blackhole 
rapidly as long as there is enough mass in the initial sterile neutrino halo. In the following, 
we assume that all sterile neutrinos fall into the blackhole and decay into active neutrinos 
and photons, and so Mbhj = Msit^) + Mbh,o ~ M^oe"'"^"*''. 



The photons emitted by the original decaying sterile neutrinos provide energy to the 
gas in the protogalaxies by Compton scattering. The optical depth of a decayed photon in 



- 6 - 



the b ulge is r = UeCTTRe, where Re is the J-band effective bulge radius (IMarconi and Hunt 
20031 ) . (Jt is Compton scattering cross section and rig is the mean number density of the gas. 
In equilibrium, the heating rate is equal to the cooling rate by Bremsstrahlu ng radiation A^, 
Recombination and adiabatic expansion A^. We have (IKatz et al.lll996l ) 



L(l 



Asontr^-' + Anon'X-' 1 + 



T 



106 K 



V + PV^'^Cs, (15) 



where A^o = 1-4 x 10"^'^ erg s"\ A/jo = 3.5 x 10^^^ erg s~^, c^, P and V are the sound 
speed, pressure and total volume of the gas within Re respectively. The Mbhj — cr relation 
can be obtained by using Eq. (15) and the Virial theorem numerically. Nevertheless, we 
first illustrate the idea by obtaining analytic relations in two different regimes. Suppose 
Ms{tb) > 10^ Mq] if r ^ 1, the resulting temperature is above 10^ K and the total cooling 
rate is dominated by Aq. For r < 1, the resulting temperature is lower and the total cooling 
rate is dominated by A^ and A^ (see Fig. 1). 



In the optically thick regime, r 1 and A^ A^ + Aji, and we get 



kT 



^ \ -2/3 



7 



1/3 



(16) 



where rrig is the mean mass of a gas particle. By using the Virial theorem kT = fiGMBrrig/SRe 
where Mb is the effective bulge mass of the protogalaxy within R^, and substituting L ^ 

MsTc^/2, we get 

2^1/2 j3/2^3/2g-r3.* /Mb 



5/2 



37/6(47r)V3rc2 [Re) ' ^^^^ 
where 7 is the adiabatic index of the gas and /i ~ 1 is a constant that depends on the 
density distribution of the protogalaxy. As time passes, the energy gained by the gas would 
decrease gradually and the mass distribution at the center would change slightly also. If 
the supermassive blackhole was formed when the galaxy formation was nearly completed 
(tb = 10^6 — 10^^ s), the ratio Mb /Re and the velocity dispersion do not change significantly. 
According to Eq. (17), the ratio Mb/ Re is fixed by Ms and F. By using the Virial theorem 
again and assuming spherical symmetry, one can relate this ratio with the final bulge velocity 
dispersion after the supermassive blackhole is formed, 

GMb 



a 



f2 



Rp 



where /2 ~ 1 is a constant that depends on the mass distribution at present. Combining 
Eq. (17) and Eq. (18), we get 



M. 



BHJ 



37/6(47^)1/3/1/2^02 



-a 



(19) 
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In the optically thin regime, r < 1 and Ar + ^ A^, the total power absorbed by 
the gas in the protogalaxies within Rf. is ~ LrieaxRe- If the cooling rate is dominated by 
Bremsstrahlung radiation, in equilibrium, we get 

kT 



By using the Virial theorem and Eq. (18), we obtain 

2ABoe-^^''*'' f h 



M, 



BHJ 



Skrrif 



1/2 



a 



(20) 



(21) 



If the cooling rate is dominated by recombination, in equilibrium and for T < 10^ K, we get 

fLarRe^''^' 



kT = k 



V Aon, 



V 



and 



2ARoe- 



BHJ 



fi 



3km. 



3/10 



a 



(22) 



(23) 



Therefore, Mbhj and a are closely related in both optically thick and thin regimes. 



We use 500 random data in the ranges of r = 0.01 — 10000, /i = 0.6 — 3, /2 = 0.6 — 3, 
= lO^s _ 1017 s, lO^M© < Mb < W^Mq and 0.1 kpc < i?e < 10 kpc to generate the 
values of Mbhj and a by using Eq. (15) an d (18) (see Fig. 2). B y using V^y = (5± 1) x 10"^^ 



which solves the cooling flow problem ( Chan and Chu 200?! ) and T^,y = 128r, we get 
'Mbhj 



log 



Mrr 







(3.98 ± 0.29) log 



200 km s- 



(8.04 ±0.20), 



(24) 



which agrees with the recent observation: a = 4.02±0.32 and f3 = 8.13±0.06 (ITremaine et al. 
2002h . 



As an order of magnitude estimate, sin 6' ~ mo/iTis, where m^) is the neutrino Dirac 
mass. In the standard s ee-saw mechanism, mp is the geom etric mean of the light neutrino 
mass-scale m,^ and (IMohapatra and Senjanovid Il980l ). Therefore, sin6' ~ -^myjms. 
From Eq. (1), we get 

.4 ^^^^ 



10" 



-23 



1 eV/ Vl keV 



For Fa. 
(m, > 10 keV). 



10 ^'^ s \ ms ~ 30 keV for rriy ~ 1 eV, which is consistent with our assumption 



^If we consider the major decaying channel is 'Sv, the decay rate obtained in this paper should 

correspond to T^^. 
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4. Discussion and summary 



We assume the existence of a degenerate neutrino halo (m^ ~ keV) at the proto-galactic 
center with a p arameter, T^^, which is universal and can be inferred from observation of 



cluster hot gas (IChan and Chul 120071 ). Without any assumptions of the protogalaxy and 
existence of any intermediate mass blackohles, a ~ 4 and /3 ~ 8 require the total decay rate 
to be = ( 5 ± 1) X 10~^^ s~\ wh ich is consistent with the observational data from cooling 
flow clusters (jChan and Chul 120071 ). Also, we have reviewed several models to explain the 
Mbhj — 0" relation. These models require several assumptions or free parameters which 
may not be true for all galaxies. For example, King's model assumes the isothermal density 



profile (p 



for all galaxies during all the time of the blackhole formation. If the density 



profile changes into the form p 



then a oc ^/r which is not a constant. This problem 



also exists in the Ballistic model which is based on the isothermal distribution of matter. On 
the other hand, in the self-similar model, there are several free parameters which are model 
dependent. Also one cannot obtain the proportionality constant f3 of the Mbhj — o^ relation. 

We have considered a wide range of r = 0.1 — 10000, /i = 0.6 — 3 and /2 = 0.6 — 3 
encompassing almost all possibilities in galaxies. In our model, we assume that the super- 
massive blackhole was formed at the epoch when the galaxy formation was nearly completed 



10 



16 



10 s) SO that the velocity dispersion does not change significantly between tf, 



and present. Therefore our result is valid only for supermassive blackholes formed nearly at 
the end of the galaxy formation, same as in the Super Eddington Accretion and Ballistic 
model. The assumed existence of a decaying sterile neutrino halo inside each galactic center 
provides enough mass to form the supermass ive blackhole. It can also solve the cooling 
f low problem in clusters fjChan and Chul 120071 ) and explain the reionization of the universe 
( iHansen and HaimanI l2004l ). all with the same decay rate T^i, = (5 ± 1) x 10~^^ s~^ and 
m. > 10 keV, which are consistent with the standard see-saw mechanism. 
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Fig. 1. — Top: The cooling rates by Bremsstrahlung radiation (dashed hne), recombina- 
tion (sohd hne) and adiabatic expansion (dotted hne) versus r in Eq. (15). Bottom: The 
temperature of the gas in a protogalaxy versus r. We used Ug = I cm~^, = 1 kpc and 
L = 3 X 10^''^ erg s^^ 




Fig. 2. — log Mbhj/Mq versus log (T200 for 500 random data, where (7200 = (7/200km s ^. We 
used r = 0.01 - 10000, /i = 0.6 - 3, = 0.6 - 3, = 10^^ - 10^^ s'^, Mb = lO'' - IO^^^q 
and Re = 0.1 — 10 kpc. The best-fit hne in the figure corresponds to slope a — 3.97 ± 0.14 
and intercept p = 8.11 ± 0.12. 



